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How to use the oven to bake delicious food is the most concerned problem of the designers and users of the oven. 


For this intent, this paper analyzed the heat distribution in the oven based on the basic operation principles and 


proceeded the data simulation of the temperature distribution on the rack section. Constructing the differential 


equation model of the temperature distribution changes in the pan when the oven works based on the heat radia- 


tion and heat transmission, based on the idea of utilizing cellular automation to simulate heat transfer process, 


used ANSYS software to proceed the numerical simulation analysis to the rectangular, round-cornered rectangular, 


elliptical and circular pans and giving out the instantaneous temperature distribution of the corresponding shapes 


of the pans. The temperature distribution of the rectangular and circular pans proves that the product gets over- 


cooked easily at the corners and edges of rectangular pans but not of a round pan. 
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Introduction 


When baking food in the oven, we all expect that the 
food looks nice, smells sweet and tastes delicious. How- 
ever, in our real life, the phenomenon of burnt food is 
quite common. In this paper, the problem to be solved is 
how to build a model to find the most appropriate pan, so 
as to make the food cannot be burnt easily. 

Firstly, this paper constructed a model to ensure the 
inner temperature distribution of the oven. The rise of 
internal temperature of the oven is mainly caused by the 
heat radiated by electric heating element. The oven has 
the impacts of the heat radiation, heat transmission, heat 
convection. 

The next, in order to better illustrate heat distribution 
on the edge and at the corner of the pan, we could cut out 
a plane to construct a two-dimensional model to consider 
the heat distribution inside one plane only in order to 
analyze the heat distribution of the whole three-dimen- 
sional model. But considering the difference of the 
shapes of the pan, range from rectangular to circle, in- 
cluding those graphics like rounded rectangle or elliptical 
which are derived from the graphics mentioned earlier, 
can make disparate influence on the heat distribution. So 
we have to do the analysis to all of them to confirm the 
optimum shape. In this process, the paper used the idea 
of cellular automation and analyzed with the help of 


computer simulation results. 


The analysis of the heating process in the 
oven and numerical simulation of the heat 
distribution on the pad section 


The three forms of thermal transmission are heat con- 
duction, heat convection and heat radiation. Due to the 
fact that the oven is made with insulating materials, there 
will be no convection, which implies heat convection 
needn’t to be considered. So the major form is heat radia- 
tion plus the heat conduction when heating products. 


The construction of the model of the heat radiation in 
the oven cavity 


Radiation heat transfer process is a process of heat 
balance in which the heat rays radiation and reflection 
transfer heat repeatedly. In this process, the final outcome 
is that the subject with higher temperature transfers heat 
energy to one with lower temperature. Therefore, the 
radiation heat transfer is related to the physical geometric 
properties of two products. 

In order to gather the temperature in oven cavity at 
any position, any time when heating, use the heat radia- 
tion equation to develop the oven-heating model. Sup- 
pose the temperature of the upper and lower heating 


pipes achieve 7; before heating, two layers of pans are 
arranged in the oven and the initial temperature is T), and 
the heat radiation of the oven is shown in Figure 1. 
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Fig.1 The heat radiation of the oven 


According to Stefan-Boltzmann Law", the equation 
for calculating the intensity of heat flow of heat radiation 
transfer is 


q, = fo(T; - Ti) (1) 

In which, q, is intensity of the radiation of heat flow, 
and the unit is W/m’. 

E= &'&, &, &stand for the degrees of blackness of 
objects 1 and 2 respectively. 

fis average angle coefficient related to the shapes and 
sizes of two objects radiate each other. 

cis Stefan-Boltzmann constant, the unit is W/(m?-K"). 

Tı is the initial temperature of the upper and lower 
heat pipes in the oven, while, T, is the initial temperature 
of the pan in the oven. 

The partial differential equation of Fourier’s law of 
heat conduction is ”!: 


V+ pe, = =0 (2) 


The heat Q = q, + qe, qc is the heat of the heat conduc- 
tion. Ignore the temperature changes on the pan surface 
therefore the qe does not exist. From the equations (1) 
and (2) we can get 

oT ô T ô T ô, T 


Pan “ata ya tae (3) 


In which, p is the density of the material, c, is the spe- 
cific heat capacity of the material, 2 is the thermal con- 
ductivity of the material. T = T(x, y, z, f) is the tempera- 
ture, which is a function related to the time ¢ and position 
(x, y, Z). According to this model, the T(x, y, z, £ inside 
the oven cavity can be obtained at any time anywhere. 


The numerical simulation of the heat distribution of 
the oven and pad section "! 


This paper has done the simulation for thermal distri- 
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bution of the rack section by the ANSYS, and the simula- 
tion results have been shown as Figure 2. 


Fig. 2 The diagram for simulation of heat distribution in the 
rack section 


The numerical simulation of the heat transmission of 
the pan and the heat distribution 


After developing the heating model inside the cavity, 
the temperature T(x, y, z, 4) can be known with the help 
of the results solved. When put pan into the oven, as the 
cavity temperature above the pan is higher than the initial 
temperature of the pan, the heat in the cavity will be 
transferred to the pan by air and then proceed the heat 
transmission above the pan. 


The process analysis of the heating of an iron wire in 
the oven 


First of all, simplify the heating process inside the 
oven at condition of one-dimensional, we can take the 
pan as an iron wire that is being heated. Assume before 
putting the baking pan into the oven, the temperature of 
the oven has been heated to a constant temperature To. 
Then in the case of one dimensional, the heat sources of 
the pan are on the endpoints. We can respectively image 
as two candles at both ends of the wire heating the iron, then 
proceeding heat conduction, which is shown in Figure 3. 
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Fig.3 Heating process ofan iron wire 
Then, use the Fourier’s heat transmission model ™! to 
analyze the heating process in one dimension. Set the 
length of the iron as L, T(x, t) stands for the temperature 
of the iron in the location x of the moment ¢, and 0 <x < 
L. So the heat conduction equation in one dimension is"! 
T, =kT, 
u(0,x)=T(x) Vx e[0,L] (4) 
T(t,0)=T(t,L)=T) Wx>0 


Heat conduction of the pan 


Analyze the heat conduction process of the pan on the 
basis of the one-dimensional heating process. Due to the 
existence of temperature diversity between upper and 
lower surface of the pan, there exists heat conduction 
between upper and lower surface, and meets the heat 
conduction equation 

O01 oT oF 
=k( + + 
ôt ax a a? 


) (5) 


in which, k= 4, . 
Pĉp 
The temperature distribution of pan on the upper and 
lower surface can be obtained based on the solution of 
thermal radiation process. For each surface, the tempera- 
ture is different in different positions. Therefore simplify 
the model to the heat conduction model from point to 
point on the top and bottom surfaces, as shown in Figure 
4. And then it can be seen as an one-dimensional heat 
conduction model !*!, 


T(0,x) =Tp (6) 
T(t,0) = T,' 
T(t) =7,' 

Then adding the temperature of all the points, we can 
get the final temperature distribution of the whole pan. In 
which To is room temperature. T,’ is the temperature of 
corresponding point on the upper surface after heat radia- 
tion. 7,’ is the temperature of corresponding point on the 
lower surface after heat radiation. The process can be 
seen in Figure 4. 


The numerical simulation of the distribution of the 
different shapes of the pan based on the cellular au- 
tomation 


The simulative basic principles of the cellular auto- 
mation 


Cellular Automation '! is an analysis model of non- 
numerical algorithm of space using simple codes and 
non-cell reproduction mechanism, and it can well simu- 
late the temperature distributions in heat transmission 
and diffusion. The structure of the cellular automation is 
shown in Figure 5. 
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Fig. 4 Thermal transmission of pan 


Cellular and its state 


Fig.5 The structure of the cellular automation 


The model analysis of the iron wire in one-dimen- 
sional 

When the model is simplified in one- dimensional, we 
can regard it abstractly as a problem, which has an iron 
wire whose one side or two sides receive heat and then 
transmit. According to Fourier’s equation of heat trans- 
mission in one dimension, 

U (t,x) = a Uy, (t >0,-0<x< +00),u(0, x) = 5 (x) (7) 

we can get the temperature of any time ¢ at any point x 
on the wire and solve it with integral transform method. 

Make the Fourier’s transform to x, u, and 6 


+0 
a(t,2)= Flu(t,x)]= f ult dexptizg jae (8) 
Suppose x > 0,u,ué — 0 , use the properties of the 


Fourier’s transform, we can get the initial problem of 

ordinary differential equation shown below 
ae =a’ l- ial = a" x? u 
dt 

u |,-9= 0 


(9) 


We are able to solve it using the methods of solving 


ordinary differential equation and get the solution 


as a(t, A) = expt a°x’t}. 


Make the Fourier inverse transform and get the basic 
solution of the Cauchy problem of the one-dimensional 
heat transmission equation 


u(t, x) = — f expt iAx\da 


+00 
-> f expt a’x t ix dA (10) 
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Then apply the cellular automation model in computer 
to simulate as shown in Figure 6. 

After the cell receives the heat from heat source radia- 
tion, it delivers it to the neighbouring cell b and b does 
the same to c, ..., as this process goes on, the heat could 
be spread on the whole iron wire. While the number of 
cells is enough, i.e, the time and space is discrete enough, 
we could gain the temperature distribution on the linear 
wire by computer simulation, as shown in Figure 7. 


Fig. 6 Cellular automata simulation process in one-dimen- 
sional case 
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Fig. 7 The temperature distribution of a wire 


The model analysis and simulation of different shapes 
of the iron plane in two dimension 


When we use the two-dimensional model to analyze 
the heat transmission on the pan, it is the same as to take 
the pan as a section so the heat spreads from the outer 
edges of the section to the center. Same as in one dimen- 
sion, we could get the temperature at any time ¢ and any 
point (x, y) on the section by Fourier’s equation of heat 
transmission"). 


1 i x+y ee 
U(t, x, y) =| — = | exp] - ————— 11 
(£x, y) (z ) 1 rer | (11) 


Once again apply the cellular automation in com- 
puter to simulate, the discrete data applying cellular au- 
tomation with different shapes of pans are shown below. 

(1) Rectangular pan 

When the pan model is rectangular, Figure 8 shows 


the simulation of the cellular automation. 

By analysis, the temperatures of a and b are the tem- 
perature inside the oven, i.e, the relative outside temper- 
ature to the pan, c and d are the temperatures on the pan. 
We could know that the product would be overcooked by 
analyzing the heat source of the cell c and d, they could 
not only gain the heat from cell a but also from b, even 
from other “neighbours”, therefore when the pan model 
is rectangular, cells in the corners have higher tempera- 
tures than others as they could be provided with heat from 
much more “neighbours”. The shape producing by the 
computer simulation through ANSYS is shown in Figure 9. 


Fig. 8 The cellular automata simulation process of two-dimen- 
sional rectangular pan 


Fig. 9 The temperature distribution of two-dimensional rec- 
tangular pan 


From the figure, we could see directly that the temper- 
ature distribution of the outer edges of the pan is uneven. 

(2) Circular pan 

While the pan model is circle, as the points on the 
edge of the circle have no difference, the heat distribution 
of the outer edge is even, as shown in Figure 10, b and 
c, as they are identical, they have same amount of the 
neighbouring cells to provide them with heat, and the 
simulation results are in Figure 11. 

(3) Elliptical pan and rounded rectangular pan 

When the pan models are elliptical and rounded rec- 
tangular, applying the same method to cellular automa- 
tion simulation and get the results of computer simulation, 
as shown in Figure 12 and Figure 13. 
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Fig. 10 The cellular automata simulation process of two-dimen- 
sional circle 


Fig. 11 The temperature distribution of two-dimensional circle 


Fig. 13 The temperature distribution of two-dimensional round- 
cornered rectangle 
From the simulation outcomes, we can see why the 


edges of food on the rectangular pan are easier to be over- 
cooked than those on circular pan. As the shapes of the 
pan evolve gradually from rectangle to circle, its temper- 
ature distribution on outer edges intends to be even. 


Why would the rectangular pan overcook 
the bread? 


Any point among the pan absorbs heat from two as- 
pects, including the heat radiation energy of cavity and 
the heat radiation energy from the outer edge of the pan. 
Suppose the heat distribution of the cavity is even, the 
temperature diversity between point A and point B is 
mainly resulted from heat radiation from the edge of the pan. 

As shown in Figure 14, take two points A and B, 
where A and B have the same distance from the edge in 
the pan, and make unit circle separately. The temperature 
of the edge of pan is identical because the heat distribu- 
tion of the cavity is even. Discretize the point of the edge 
of pan, and regard these points as the source of thermal 
transmission, thermal radiation of any point in the edge is 
identical. From the geometry we can know /; + h > J, so 
the quantity of the source of the heat gathered around 
point A is bigger than B, which implies the thermal radia- 
tion of point A gathered from the pan is bigger than B. 
Hence the temperature in the edge of the pan is higher 
than others. That’s why the product in the outer edge of 
the pan is easy to be overcooked. 


Fig. 14 The sketch map explaining the high temperature in the 
four corners 


Conclusions 


From the heat distribution model and numerical simu- 
lation results in this paper we know that, due to uneven 
heat distribution at the edges and corners of rectangle pan, 
it would be easy to burn the food. However, if using cir- 
cular pan, as a result of even heat distribution on the en- 
tire pan, the food would not be burnt easily. But given the 
truth that the ovens used in the actual life are mostly rec- 
tangular, using rectangular pan can make full use of the 
oven space. We can bake more food in one time. On the 
contrary, circular pan cannot make good use of the oven 


space. 
In actual use, we can choose rounded rectangular pan. 
It not only can ensure that food would not be burnt easily 
at the edges and corners of the pan, but also can make full 
use of oven space. Of course, the specific choice of the 
shape of the pan depends on the user’s preferences. If you 
just want to bake delicious food, choose the circular pan! 
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